Magnetic bit stability: Competition between domain-wall and monodomain switching 
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We numerically study the thermal stability properties of computer memory storage realized by 
a magnetic ellipse. In the case of practical magnetic random-access memory devices, the bit can 
form a spin texture during switching events. To study the energy barrier for thermally-induced 
switching, we develop a variational procedure to force the bit to traverse a smooth path through 
configuration space between the points of stability. We identify textured configurations realizing 
domain- wall propagation, which may have an energy barrier less than that of the corresponding 
monodomain model. We contrast the emergence of such micromagnetic effects in thermal versus 
field-induced switching. 



The mechanisms of tunnel magnetoresistance [1 and 
spin-transfer torque [2 offer a simple and localized 
method in which to read and write computer memory by 
applying electric currents rather than magnetic fields [3] . 
The corresponding devices have the potential to improve 
write energy and speed and decrease bit area. Equally as 
important is the lifetime of the bit; that is, the robustness 
of the bit against thermal fluctuations between the "on" 
and "off' states. This is furnished by the energy barrier, 
£^5, separating the points of bimodal stability and typi- 
cally characterized by the dimensionless thermal-stability 
parameter, A = ^^^//c^T, wherein /c^ is Boltzmann con- 
stant and T ambient temperature. 

There are principally two classes of methods to experi- 
mentally extract this thermal barrier: stochastic and field 
induced. In stochastic methods, the thermal barrier is ex- 
trapolated from the thermal fluctuations of the magnet 
between the bistable points. Because these fluctuations 
are of the order of the bit lifetime at room temperature, 
switching is in practice assisted by a large temperature, 
a spin- transfer torque, or a magnetic field [4 . However, 
at large temperatures, it is becoming difficult to extrapo- 
late the system parameters, such as magnetic saturation 
and anisotropics, making this method unreliable. Like- 
wise, although the thermal stability can be extrapolated 
faithfully at low voltages upon the application of a spin- 
transfer torque, at high voltages the local temperature 
may be too high to rely on the equilibrium parameters of 
our system due to the injection of "hot" electrons and the 
associated Ohmic losses. For deterministic field-induced 
switching, one applies an external magnetic field in the 
easy plane (Sj . By determining the magnitude of the mag- 
netic field at which the system is no longer bimodally sta- 
ble, the energy barrier can be ascertained. Following this 
procedure for all directions in the easy plane, one obtains 
a Stoner-Wohlfarth astroid curve [6 . In order to be vi- 
able, all these methods generally assume a single-domain 
model of the magnet. In this model, there is only one pla- 
nar path through configuration space between one point 
of stability and the other. However, when the character- 
istic length scale of the spin texture becomes comparable 
to the major axis of the elliptical bit, a richer configura- 
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FIG. 1. Application of the local pinning fields (thick arrows) 
and schematic of the resulting magnetic textures (thin ar- 
rows): (a) Center pinning, (b) left off-center, (c) simultaneous 
left-off-center and right-off-center, and (d) simultaneous cen- 
ter and left-off-center. We find that (b) produces the lowest 
energy barrier. 



tion space becomes energetically accessible, opening the 
possibility for a textured path with a lower energy barrier 
between the stable points. As a result, an extrapolation 
of the thermal barrier can be ambiguous. 

In our analysis, we consider a magnetic ellipse of typi- 
cal dimensions and material properties as those expected 
to be used as memory. For these values, we find that 
the exchange energy is of the order of the demagneti- 
zation energy. By applying slowly rotating local mag- 
netic fields, which guide some variationally preselected 
points along the major axis (see Fig. [T]), we force the 
bit into a series of textured configurations through phase 
space that continuously evolve one stable minimum into 
the other. We find a path between the bistable points, 
resembling domain- wall propagation, whose energy bar- 
rier is smaller than that of a single domain. For a con- 
crete numerical demonstration, we model our bit as a 
V — 150 X 50 X 2 nm^ thin-film ellipse with saturation 
magnetization of Ms = 1100 emu/cm^ and exchange 
stiffness A = 2 /ierg/cm. The crystalline anisotropy is 
neglected, and we consider only this single layer. Our 
axes are oriented such that the xy plane is in the plane 
of the ellipse and the easy (major) axis is along x. See 
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FIG. 2. Energy as a function of angle for different pinning 
prescriptions: (a), (b), (c), (d) of Fig. [l] monodomain (e), 
and left off-center pinning while initializing the bit at ^ = 180° 
(f). Insets (i)-(iii) show rux at several critical points along the 
curves, with red corresponding to positive and blue negative 
magnetizations. 

Fig. [l] The energy of the bit is written as sum of the 
exchange and demagnetization contributions: 

F = E^c + ^dm , (1) 

where E^c = ^/ dV{\/m)'^, m is the unit magnetiza- 
tion vector, and £^dm is the dipole-dipole interaction. Be- 
cause the applied magnetic field is used only to simulate 
a texture through phase space, in what follows we ig- 
nore the energetic contribution of the applied field. In 
the single-domain limit, E^c = and = {^x^x + 
NyiJiy -\- Nzm1)VM^ ^ where the demagnetization factors 
Nx + Ny -\- Nz = 27r. The energy barrier is then simply 
given by £^5 = {Ny — Nx)VMg^ corresponding to the xy- 
plane magnetization rotation. We define the character- 
istic length of the texture to be ^ = V^/i^y ^x)M^. 
Because for our structure £ ^ 30 nm, less than the ma- 
jor axis of our ellipse, micromagnetic (i.e., spin-texture) 
effects can be important. The latter are modeled using 
the LLG Micromagnetics Simulator j7j- 

To find a continuous path through configuration space 
from one stable point to the other, we apply a single 
or a pair of localized pinning fields that force the mag- 
net into a particular configuration, as sketched in Fig. [T] 
Initially oriented along the easy axis, we slowly rotate 
the pinning field in the xy plane, allowing the (twisted) 
texture to equilibrate, and calculate the internal energy. 
In order to variationally find an optimal path, we stud- 
ied four different methods of applying the pinning field, 
which are depicted in Fig. [l] Calculating the internal 
energy as a function of angle of the pinning fields, we 
are able to plot the energy of the system as the bit flips. 
The principal results of our work can be found in Fig. |2j 



Under application of a single off-center pinning, the en- 
ergy barrier, A ^ 90, is less than that of a monodomain, 
A ^ 110. Since the thermal energy barrier is lowered by 
micromagnetic effects, the bit will follow a textured path 
in phase space when thermally flipping. Underlying the 
seemingly small change of < 20% between energy barri- 
ers is thus a fundamental change in the mode of thermal 
switching. The Arrhenius thermal transition rate oc 
is, furthermore, greatly reduced. Upon examining the 
texture dynamics, one finds that the left off-center pin- 
ning field nucleates a domain wall, while rotating the 
pinning field pushes the domain wall to the right along 
the X axis. This is most apparent at about 145°, where, 
as illustrated in inset (i) of Fig. |2j the two domains each 
occupy approximately half the area of the magnet. At 
this point, the total internal energy is at a maximum, 
held in place only by the pinning field. After moving the 
domain wall past this point, there is a precipitous drop 
in internal energy as result of the domain-wall propa- 
gating to the right side of the magnet. This rapid drop 
in energy reflects a hysteretic process, wherein initially 
orienting the magnet at 0° or 180° with respect to the 
X axis, the energy will respectively follow the (b) or (f) 
curve of Fig. |2j We conclude that it is the domain-wall 
nucleation and propagation process that thermally flips 
the bit in our example. 

To investigate the micromagnetic effects in the field- 
induced extraction of the thermal barrier, we apply a 
uniform magnetic field antiparallel to the direction of 
initialization of the bit. Our magnetic bit flips at be- 
tween 600 and 650 Oe, which is close to the value pre- 
dicted by the single-domain anisotropy: 680 Oe. Closely 
examining the texture, we indeed see the magnet flips 
like a monodomain. Analogous to the case of thin rect- 
angular magnetic strips, wherein the coercivity field in- 
creases as the lateral width decreases, we suggest that 
the energetic bottleneck for domain-wall switching un- 
der application of an external field is at the tip of the 
bit where the domain wall nucleates. In contrast, the 
Arrhenius barrier for thermal domain- wall switching cor- 
responds to the domain wall at the center of the mag- 
net. We can thus distinguish three principal scenarios 
subject to the characteristics of the bit: First, the bit 
flips like a monodomain under both thermal fluctuations 
and application of a coercivity field (as in the case of 
infinite stiffness), illustrated schematically in Figs, [sj^c) 
and (d), respectively. Second, the bit flips thermally 
via domain- wall propagation [Fig. |3|b)] and like a mon- 
odomain under application of a coercive field [Fig. |3|d)] 
(as in our example). Third, the bit flips through domain- 
wall propagation both thermally and under application 
of a coercive field [Figs. [3ja),(b)]. We find this in our 
simulations on more elongated ellipses, e.g., with dimen- 
sions 400 X 80 X 2 nm^ (and the same material param- 
eters). The stochastic or field-induced characterization 
of a physical bit could correspondingly be significantly 
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FIG. 3. Schematic illustration of the "bottleneck" magnetic 
configurations for switching under the application of magnetic 
field, (a) and (c), and during thermal fluctuations, (b) and 

(d). 



affected depending on the appropriate regime. 

Since the exchange energy scales linearly with bit 
thickness while dipolar energy quadratically, upon in- 
creasing (decreasing) the thickness the texture effects are 
enhanced (reduced), and, therefore, the disparity in the 
thermal barrier between the single-domain path and the 
domain- wall should increase (decrease). Particularly we 
take a thickness of 3 nm (1 nm) and find that the mon- 
odomain is > 20% (< 1%) larger than the micromagnetic 
thermal barrier. 

We may expect random impurities and edge roughness 
to accentuate texture of the device and therefore affect 
the relative size of the thermal barrier in micromagnetic 
versus monodomain models. As an example, we model 
roughness as a crude edge modulation with wavelength 
20 nm and amplitude 10 nm, as depicted in the left half 
of Fig.|4] While this is perhaps larger than the character- 
istic edge roughness of a typical device, we take this as a 
somewhat exaggerated case. Such edge modulation may, 
furthermore, be used to represent a generic structural 
or material disorder. We find that the thermal energy 
associated with both modes (thermal and field-induced) 
of switching is diminished by disorder, while the rela- 
tive difference between the single-domain and textured 
behavior, < 20%, is essentially unaffected [see Fig. [4ja)]. 
Further distorting the edge (right half of Fig.|4|, the ther- 
mal barriers become nearly equal [see Fig. [4]^b)]. By way 
of roughening the edge, the aspect ratio of the ellipti- 
cal bit has decreased, decreasing the xy demagnetization 
factors and thus favoring monodomain switching. A sys- 
tematic analysis of geometric and material imperfections 
falls outside the scope of this work. 

As a result of our variational optimization of the path 
taken through the configuration space, there may be a 
trajectory of states with even a lower associated ther- 
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FIG. 4. Energy profile of bit with (a) edge roughness of ampli- 
tude 10 nm and wavelength 20 nm and (b) further enhanced 
roughness, where the dashed curve shows single-domain and 
the solid curve micromagnetic switching obtained under ap- 
plication of a local pinning field to the left off-center. The 
corresponding equilibrium bit textures are shown below. 



mal barrier. In other words, the domain-wall type of 
propagation discussed here produces the upper bound 
for the thermal barrier. Finding at least one micromag- 
netic path, however, whose thermal barrier is smaller 
than that of the monodomain, suggests a qualitatively 
different scenario of thermal flipping and a lower thermal 
barrier, compared to the monodomain. We find such ef- 
fects, with regard to both thermal stability and astroid 
diagrams, ubiquitous in practical devices. 
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